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It is pointed, that effects of refraction of electromagnetic radiation in the medium, 
formed by the magnetized vacuum, become essential already for relatively soft 
photons, not hard enough to create an electron-positron pair, including those be- 
longing to soft gamma-. X-ray, optic and radio- range, if the magnetic field B 
exceeds the critical value of Bcr = rn? je = 4.410^3 Gauss. 



Three leading terms in the asymptotic expansion of the one-loop polarization op- 
erator in a constant magnetic field are found for B > > Bcr, and the corresponding 
refraction index is shown to depend only on the propagation direction of the pho- 
ton relative to the external field. It is established, that the refraction index for one 
of polarization modes unlimitedly grows with the field, while the other is saturated 
at a moderate level. The photon capture effect is extended to soft photons. 



The results may be essential in studying reflection, refraction and splitting of X- 
rays, light and radio waves by magnetic fields of magnetars, as well as in considering 
emission of such waves by charged particles . 
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1 Introduction 



Although it is long since the refracting and birefringing properties of a strong 
magnetic field in the vacuum have been realized, up to now their only essen- 
tial consequences considered in a realistic astrophysical c ontext remai n the 
photon splitting effect ^ and the effect of photon capture l^l^|4 | 5 | 6 [7|8| j^^ie 
both effects are currently discussed mostly as applied to electro-magnetic ra- 
diation belonging to the 7-range. They depend crucially on the deviation of 
the photon dispersion curve from the customary shape it has in the empty 
vacuum, fcg — |kp, where ko is the photon energy, and k is its momentum. 
As long as one restricts oneself to considering the magnetic fields B below 
the Schwinger critical value B < Bcr = [m? /e) = 4.4 • 10^^ Gauss, where m 
and e are the electron mass and charge, resp., the only essential source of this 
deviation is the singular behavior of the polarization operator Tl^^{k) near the 
creation thresholds of mutually independent electron and positron on Landau 
levels n, n' by a photon (the cyclotron resonance) or a still stronger sin- 
gular behavior of II^j, near t he poin ts of a mutually bound e+e^-pair (the 
positronium atom) formation To reach (at least the lower of) these 

positionss the photon should belong to the 7-ray range, with its energy above 
or of the order of 1 Mev. For this reason the effect of photon capture, with its 
transformation into electron-positron pair, derived from the singular behavior 
of n^i,(/c), concerns mostly the 7-quanta, as long as their propagation in a 
pulsar magnetosphere of traditional pulsars is concerned. It was estimated, 
that the fields about B — O.IBcr are sufficient to provide this effeciP^and to 
protect the positronium atom, into which the captured 7-quantum is trans- 
formed, against ionization by the acc elerating electric field in the polar gap 
and by the thermal photons EEEEEl. 

Also the Adler effeclPof photon splitting 7 ^ 77 for such fields is usually 
discussed relating to 7- quanta . There are two reasons why, again, 

the 7-range is important. The first one is that the photon splitting becomes 
possible in the magnetic field, because the deviation of the dispersion curve 
from the fcg = law opens a kinematical aperture for this process - the wider, 
the stronger the deviation ( and the deviation is strong near the thresholds). 
Besides, there is a strong birefringence for the photons in the 7-range, since 
only one eigenvalue K2{k) of the tensor II^i, is singular near the lowest (n = 
n' = 0) threshold, while the other two eigenvalues ki,3(A:) remain finite, until 
the next thresholds (n = 0, n' = I or n = 1, n' = 0) are reached. This 
implies that the photons of only one polarization mode are essentially affected 
by the medium. This birefringence leads to polarization selection rules in 
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the photon spUtting process, which are well pronounced. The second reason 
is dynamical. The matrix elements of the photon splitting are subject to 
the same resonant behavior near the thresholds as the polarization operator. 
The aforesaid explains why mainly the 7-rangc is first to be affected by the 
magnetized vacuum. 

The situation changes considerably after one passes to superSchwinger 
magnetic fields B Bcr, expected to be existing in soft gamma-ray repeaters 
and anomalous X-ray pulsars, whose magnetic field is estimated that large ( 
I refer to the talk, given by Prof. A.Reisenegger in this Workshop 'WJ also 
to El). The point is that in such asymptotic regime a linearly growing term. 




the polarization operator, thus providing an extra ( apart from the cyclotron 
resonance) large contribution into the refraction of the vacuum. 

In Section 3 we study the consequences of this phenomenon for the pho- 
ton propagation, basing on the first three leading contributions into asymp- 
totic expansion of the polarization operator eigenvalues for large B, obtained 
within the one-loop approximation. One of these consequences is a frequency- 
independent, but direction-sensitive, large refraction index for propagation 
nonparallel to the magnetic field in one ( out of three) polarization modes 
in the kinematical domain far from the threshold. The corresponding strong 
polarization- and direction-sensitive refraction occurs for electromagnetic ra- 
diation of any frequency range, including X-ray, optic and radio range. 

This study is preceded by Section 2, where exact results concerning 
the electromagnetic radiation propagation in the magnetized vacuum are de- 
scribed. These follow exceptionally from the general properties of relativistic, 
gauge, charge invariance'^-' and Onsager theorem The results, presented 
in Section 2, are valid irrespective of any approximation and the field strength, 
unless the opposite is exphcitly indicated. 

In Appendix the asymptotic expansion, used in Section 3, is derived. 

2 Exact Facts about Electromagnetic Eigenmodes in an 
External Magnetic Field 

There are three propagating eigenmodes, corresponding to the vacuum exci- 
tations with photon quantum numbers in an external magnetic field B. The 
dispersion law, i.e. the dependence of the energy kg of the quantum (or the fre- 
quency in the wave) on its momentum k, is given for each mode by a solution 
of the equation 





proportional 



one of the eigenvalues, K2, of 



(1) 
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Here k\\ and fcj^ are the momentum components along and perpendicular to 
the magnetic field B, resp., is the photon four-momentum squared, fc^ = 
+ ^1 ~ ^0- "The k's in the right-hand sides in Eqs.l^) are eigenvalues of the 

polarization operator! -*- '^ 1 "^ 1 ^ 1 

It is a general consequence of the relativistic covariance that the eigen- 
values depend upon the two combinations of momentum, specified in 
This implies, that solutions of the dispersion equations Q have the following 
general structure 

fco = fc{|+/.(fci), i = 1,2,3, (2) 

and that the direction of the group velocity in each mode v = dk^) j 9k does 
not coincide (as long as k^ ^ 0) with that of the phase velocity k/fco. To see 
this, calculate the components of the group velocity across and along the 
magnetic field B, resp., on solutions (0) of each dispersion equation Q 



-1 dnj 

dko k± dk^ k± ^ k± dfi{kj^) 

dk±_ ko dk\ ko 1 + g^^^f^j^p dfc^^ 

_ dko h 



(3) 



It follows from Q that the angle 9 that the direction v of the electromag- 
netic energy propagation makes with the external magnetic field satisfies the 
relation 

v± n dKi \ I dKi 




tane = ( 1 - ^ ) I 1 + tani?, (4) 

where i? is the angle between the photon momentum (phase velocity) and 
the external field, tani^ = k±/k^^. The following statement takes place: if 
the phase velocity exceeds the velocity of light c, i.e. if fc^ -I- fcy > feg, (or 
fiik^) < k\ in ^ ), but the group velocity Q does not, v\ + v'^^ < \ (or 
d^/i(fci)/(dfcl)^ < 0, then tan 6* < taniJ. The conditions of this statement are 
fulfilled for the dispersions laws found within the approximation-dependent 
calculations of the k's. For the superSchwinger fields, treated within one-loop 
approximation, this fact follows explicitly from equations of Section 3 below. 
Therefore, the photon tends to deviate closer to the magnetic field line, as 
compared to its momentum direction. 

It follows from the gauge invariance that 

K,(0,0) = 0, i = 1,2,3. (5) 
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This property implies that for each mode there always exists a dispersion 
curve with /i(0) = 0, which passes through the origin in the {k^ — A:|,/c^)- 
plane. However, out of these three solutions only two may simultaneously 
correspond to physical massless particles - the photons. The third one is a 
nonphysical degree of freedom, characteristic of gauge theories: in a mag- 
netic field a photon has two degrees of freedom, the same as in the empty 
vacuum. Which of the modes becomes nonphysical, depends upon direction 
of propagation and specific form of the function fi{k\) in (|2Jl. We shall fur- 
ther discuss this matter for the superSchwinger field limit in the next Section. 
Massive branches of solutions of with /i(0) > 0, may also exist, despite 
(O. For them the number of physical degrees of freedom is three, so all the 
three equations may have ph ysical solutions simultaneously (see, e.g. the 
positronium branches found in 

The refraction index Ui in Mode i is 



Unlike k^, the refraction index Ui is not a Lorentz scalar and may depend on 
two energy-momentum variables, after it is reduced onto the dispersion law 
(0). The gauge invariance property © implies that the refraction index jnj 
for parallel propagation, k^ — 0, be exactly equal to unity for the massless, 
/i(0) = 0, branches in every mode 

n| = 1. (7) 

The electromagnetic wave propagating strictly along the external constant 
and homogeneous magnetic field does this with the velocity of light in the 
vacuum c, the phase and group velocities coinciding in this case. 

If, within a certain approximation, the eigenvalue linear function 

of its arguments with the condition (0) met, the refraction index © for the 
corresponding dispersion law depends on a single combination of the photon 
energy and momentum, which is the direction of propagation d. This happens 
in a nonresonant situation, for instance, as described in the next section. 

The pol ariza tions of the modes are described in an approximation- inde- 
pendent the following relations, where e^'^ and h^'' are the electric 
and magnetic fields in the wave belonging to Mode number i — 1,2,3 

e(^)^-^fco, hW = (^xk|,), (8) 

k\ k\ " 



e 



k„ 
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e(3) = -i.of^x^V h,(-^)=-^fc,|, h||(3) = llfc,. (10) 
\k± fc|| J /c^ " " /c II 

Here the cross stands for the vector product, and the boldfaced letters with 
subscripts || and ± denote vectors along the directions, parallel and perpen- 
dicular to the external magnetic field, resp. The electric field e in the wave 
of Mode 1 is parallel to kj^, in Mode 2 it lies in the plane containing the 
vectors k, B, in Mode 3 it is orthogonal to this plane, i.e. Mode 3 is always 
transversely polarized. 

Note, that the normalization is different for each line of (|10|1 . so one may 
judge about vanishing of some components as compared to others within one 
line, but not between different lines. 

Concerning the direction of propagation, two cases are essentially differ- 
ent. If k± — 0, we speak about longitudinal propagation. If not, a Lorentz 
boost along the external (constant and homogeneous) magnetic field exists, 
which does not change the value of the latter and does not introduce an extra 
electric field, but nullifies Hence, the general case of nonparallel propa- 
gation k± ^ 0, A:|| 7^ is reduced to purely transversal propagation, k^^ = 
(in the corresponding reference frame). One should keep in mind, however, 
that the above transformation changes the photon energy kg and should be 
treated with precautions when one considers a field with curved lines of force. 

For transversal propagation, k _L B, (fcy = 0), Modes 2, 3 are transversely 
polarized e^^^'^^ _L k in two mutually orthogonal planes, e^^^ ± e'^^ while 
Mode 1 is longitudinally polarized, e*^^) || k with no magnetic field in it, 
h^^) = 0. It is expected not to correspond to a photon (depending on the 
dispersion law). 

On the contrary, for longitudinal propagation, k || B, {k± — 0), Modes 
1, 3 are transversely polarized, e'^^''^^ _L B, and their electric field vectors lie 
in mutually orthogonal planes, e^^^ _L e*^'^^ as they always do, while Mode 2 
is longitudinally polarized, e*^^^ || B, and does not contain a magnetic field, 
h^^^ = 0. This time Mode 2 is expected not to correspond to a photon, whereas 
Mode 1 is a physical electromagnetic wave, which matches the electromagnetic 
wave of Mode 3: together they may form a circularly polarized transversal 
wave because of the degeneracy property 

K^iikl-kf),0) = n,{{kl-kl),0). (11) 

The latter reflects the cylindrical symmetry of the problem of a photon prop- 
agating along the external magnetic field. 

Another remark of almost general character is in order. One might ex- 
pect the possibility of the Cherenkov radiation by a charged particle moving 
in an optically dense medium formed by the magnetized vacuum. This effect 
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- with the Cherenkov photons softer than fco = 2m - does not take place in 
known situations, however. Consider an emission of a photon by an electron 
in a magnetic field, not accompanied by the change of its Landau quantum 
number, n — n' (otherwise, that would be the cyclotron, and not Cherenkov, 
radiation). According to t he kinematical analysi s of the energy and momen- 
tum conservation, done i iJ^, (also to the studjl^ of analyticity regions of 
the one-loop photon polari zatio n operator in electron-positron plasma in a 
magnetic field, calculated i the Cherenkov photon with fco < 2m may 
only belong to the right lower sector 



in the (fcg — fcy , fc^)-plane. The substantial reason for this circumstance is 
the degeneration of the electron energy with respect to the center-of-orbit 
position in the transversal plane. No dynamical calculations, hitherto known, 
provide a penetration of photon dispersion curves into this sector. The only 
exception is a nonphysical situation due to exponentially strong external fields, 
to be mentioned in Subsection 3.2 below. We conclude that under standard 
conditions no Cherenkov emission of a photon softer than 2m is possible. 

3 Photon Dispersion in SuperSchwinger Magnetic Field 

3.1 Asymptotic expansion of polarization tensor eigenvalues 

In the asymptotic region of supercritical magnetic fields B ^ B^r and re- 
stricted energy of longitudinal motion fcg — fc^ <C {B / Bcr)m'^ , the three eigen- 
values Ki,2,3(^) of the polarization operator - provided the latter is calculated 
within the one-loop approximation as done in ' ' ' - have the following be- 
havior, derived from equations of Ref.l^ (see Appendix) 



fcg - fcjj < 0, fci > 



(12) 




(13) 



K2(fc^ -fcf,fci) 




(14) 



Ksi^Q fell , fc^) 
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Here a = 1/137 is the fine structure constant, C — 0.577 is the Euler 
constant. Eqs. p3|l . p5|) are accurate up to terms, decreasing with B hke 
{Bcr/B)\n{B/Bcr) and faster. Eq. p4|) is accurate up to terms, logarithmi- 
cahy growing with B. In K1.3 we took also the limit fc^ ^ {B / Bcr)'m^ , which 
is not the case for K2, wherein the factor exp {—k^Bcr/'2m^B) is kept 
different from unity, because it is important near the cyclotron resonance, 
as it will be explained in Subsection 3.2 below. The integral in (|14|) can be 
readily calculated, but we shall not need its explicit form here. 

The parts growingwith B in ki,2,3 we re writ ten in^l their derivation 
from equations of Ref. El is traced in detail inl^EH. The linearly growing term 
in Eq. p4|l in a different way was obtained in ^1 using a two-dimensional (one 
time, one space) diagrammatic technique developed to serve the asymptotic 
magnetic field regime. The logarithmic terms in the expressions above do not 
dominate over the constant terms unless one would like to include exponen- 
tially large magnetic fields into consideration. (That would be unreasonable 
not only because such fields are hardly expected to exist in nature, but mainly 
because their consideration is beyond the scope of Quantum Electrodynam- 
ics: the logarithmically growing terms in p3|l . (|15|l are associated with the 
absence of asymptotic freedom in QED (c/. analogous asymptotic behavior 
l^in the Euler-Heisenberg effective Lagrangian).) The derivation of all terms 
in Eqs. (|13|l . p4|l . H15|l . including those that do not grow with B, is given 
in Appe ndix us ing a straightforward method differe nt f rom the one applied 
earlier in l ^l l ^^ l The asymptotic expressions used in ^1 do not coincide with 
ours, except for the linear-in-B term. 

The limiting expressions (|13|) . H14|l . (|15|) do satisfy the exact properties 
(|TT|l and ©. 

In this paper we shall only deal with the transparency region, fcp — fc^ < 
4m^ (i.e. with the kinematical domain, where aci^2,3 are real), because we shall 
be interested in photons with fco < 2m, even fco ^ 2m, which never reach the 
free pair creation threshold fcg — fc^ = 4m^. The eigenvalue K2 (I14II has a 
singular branching point in the complex plane of the variable (fe§ — fc^) near 
the lowest pair creation threshold (fcg — k'^^)thr — ^rn? ■ Thresholds of creation 
of e"'"e~-pairs with the electron and the positron on excited Landau levels 
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n,n' ^ Q 



B 



,_B_ 

Brr 



(16) 



are shifted in the asymptotic regime to infinitely remote region. For this reason 
the eigenvalues K1.3, which are responsible for photons of such polarizations 
that can only create e+e_- pairs with at least one charged particle in an 
excited Landau state, do not contain imaginary parts or singular branching 
points in this regime. On the other hand, the eigenvalue K2 has only one 
singular branching point, corresponding to possibility of creation of electron 
and positron in the lowest Landau states by the photon, polarized as in Mode 
2. The singular threshold behavior of (|14|l near the point fcp — fcjj = 4m^ — e, 
e > 0, e ^ is 



2aBm^ 
H2[k) ~ — exp 

JDrr 



■)( 



Am^ - kl 



(17) 



As might be expected, this is the same as the behavior near this threshold 
of the exact one- loop expression for K2{k) , before the limiting transition to 
large fields has been performed. 



3.2 Propagation of eigenmodes in the superSchwinger field limit 

If Eg. ljlSI) for ki is taken as the r.-h. side of Eq. the latter has only 
one solution, which is the trivial dispersion law k^ — 0. With the relation 
fc^ = satisfied, however, the 4-potential corresponding to the electromagnetic 
field of Mode 1 becomes proportional to the photon 4- momentum vector ku, 
unless fc|| = (see 1 *^ 1 ^ -"- l ^^ l). Therefore, for nonparallel propagation. Mode 1 
corresponds only the gauge degree of freedom we discussed in Section 2, with 
no real electromagnetic field associated with it. 

Solutions of equation for the second mode i — 2 with Eq. H14I) taken for 
K2 are plotted in the Figure for three values of the field B using MATHCAD 
code. These are dominated by the cyclotron resonance (|17|) . that causes the 
strong deviation of the dispersion curves in Fig.l from the shape k'^ — { the 
light cone ). Once fc^ — > 00 near the threshold on the dispersion curves, the 
quantity k^Bcr/rn^B should be kept different from zero even in the large 
field limit under consideration. 

The behavior of the dispersion curves of Mode 2 near the threshold for 
superSchwinger magnetic fields B Bcr is the same as for the "moderate" 
fields B < Bcr, and therefore it also presents the photon capture effect for 
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Figure 1. A family of dispersion curves for Mode 2 - solutions of equation with 
Ea. ll4i taken for the r.-h. side, - below the threshold fcg — fc| = 4m^. The values 
of the external magnetic field, corresponding to the curves, are ( from left to right) 
B = 10-Bcr, 100-Bcr, 1000-Bcr- The straight line is the light cone - dispersion curve for 
B = 0. The dashed horizontal line marks the maximum to which the photon with the 
energy kg may proceed , provided that fcg < 2m. The variables along the axes are plotted 
in units of im^ IMev'^ 

photons harder than 2m, known for such fieldfl if one calculates l@J near the 
threshold fcg — A;| = 4m^ using Eq. ((T7|) as K2 to get 

tan0 = |^#^(4™^-fco' + fcf), (18) 
K|| Bm^ I II 

one concludes that the angle 6 between the external magnetic field and the 
direction of the wave packet propagation in Mode 2 tends to zero, the faster 
the stronger the field. If the photon energy ko is just a little less than 2m, 
the photon may be close to the threshold when its fc|| disappears. In this 
upper point the wave packet stops, since the group velocity length + = 
w^(l + tan^ 9), equal to fc|/fco according to the second line of Q and lfTH|l . 
disappears together with fc|| . 

As applied to the conventional pattern of a pulsar magnetosphere, already 
described in the talk by Prof. Qiac"^"^, this effect acts as following^. A curva- 
ture 7-quantum emitted tangentially to the magnetic line of force, i.e. placed 
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initially in the origin of Fig.l, then evolves along its dispersion curve as it 
propagates in the dipole magnetic field with its line of force curved, since the 
components fcy and k± are changing. The maximum value of the ordinate 
fcp — fc^ occurs at fc|| = 0, and it is the photon energy squared, k^. If the lat- 
ter is greater than 4™^, the photon may achieve the horizontal asymptote in 
Fig.l. Here its group velocity across the magnetic field, dfco/dfc_L, disappears, 
dko/dk± 0, and hence it propagates along the magnetic field and does 
not cross the threshold, since the other branch of the dispersion curve, which 
passes above the threshold, is separated by a gap from the initial branch. As 
a matter of fact, a mixed state: photon-pair is forme(3, analogous to polari- 
ton known in condensed matter physics. The massless part of its spectrum is 
presented by the dispersion curves in Fig.l. The photon gradually turns into 
the — e^-pair and exists mostly in that form when finally it is propagating 
along the magnetic lines of force. This capturing effect is important for the 
formation of radiation of pulsars with the fields B > O.lBcr, since it prevents 
- once the bindin g of the e lectron-positron pair into a positronium atom is 
taken into accounl^EEEEl _ the screening of the accelerating electric field in 
the polar gap. It may be essential for magnetars with their fields ^ 10^^ — 10^^ 
Gauss, as well. 

The new features, introduced by superSchwinger fields are that the dis- 
persion curves for Mode 2, presented in Fig.l, step aside from the light cone 
already far from the resonance region. This means that although the photons 
softer than 2m = IMev cannot proceed to the values of the ordinate in Fig.l 
higher than their energy squared (corresponding to fc|| = 0), they can still 
reach the region where the transversal group velocity dfco /d-k± becomes much 
less than unity and are thus captured to the trajectory almost parallel to the 
magnetic field. This is how the capture effect extends to the photon energies 
below the border fcg = 2rn. The cyclotron singularity at the pair-creation 
threshold in such fields is so strong that even low energy photons that are 
unable to create a pair are sensitive to it, provided that they belong to Mode 
2! 

Besides the extension of the photon capture effect to softer photons, the 
inclusion of superSchwinger fields into consideration has another impact. It 
leads to large direction-dependent refraction of Mode 2 electromagnetic waves 
of low frequency. To see this, consider the limit 

fe^ - fejj < (19) 
in Eq. H14() , which reduces to neglecting kg — k?. in the integrand in (|14|l . Then 



newhepp: Workshop SMFNS/ICIMAF - February 1, 2008 11 



(|14|l becomes 



M^) = ^«-'=i)£exp(-i^). (20) 



The exponential factor in (|20l) cannot be essential within the region p9|l . 
The dispersion equation for Mode 2 (j = 2) then has the solutions ex- 
pressing the photon energy fcg as a function of its transversal and longitudinal 
momentum: 

Ea. (|21|l presents analytically the straight line parts of the dispersion curves 
in Fig.l adjacent to the origin for various values of B. The components of the 
group velocity Q f_L,|| calculated from (|21|l are 



k± f a B \ ^ fell 



The modulus of the group velocity squared is now 

1 S ^Qg2^ 

2,2 ^ , Sir Bcr ^'oo^ 

V, +v<, = 5 — 5 , (23) 

where d is the angle between the photon momentum and the field, tan-d = 
k_L/k\\. Eq.(|23) has the maximum value of unity for the parallel propagation, 
I? = 0, in accordance with the general statement of Section 2, and is minimum 
for perpendicular propagation, ■& = 7r/2. 

The expression Q for the angle 9 between the direction of the electro- 
magnetic energy propagation and the external magnetic field in the super- 
Schwinger limit for Mode 2 becomes 

— =tan0==-iH =tan?9H . 24) 

v\\ fc|| \ Stt Bcr J \ Svr Bcr J 

Once tan0 < taniJ, in the problem, formulated above, the photon emitted 
tangcntially to curved lines of force will bend towards these lines. This relates 
to low-frequency radiation, as well. 

The refraction index © in Mode 2 for — kj <s: 4m^, B > Bcr is 

^J^B^\ ( ) 
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The refraction index obtained depends upon the direction of the photon mo- 
mentum, characterized by the angle z?, but does not depend upon its energy. 
In other words, there is no frequency dispersion in a wide range from slow 
radio waves and up to soft 7-rays with fco ^ 2m. This is a consequence of the 
fact, that in K2 only linear parts in momenta squared were as a matter of fact 
left ( correspondingly, ./(fc^) in ^ is proportional to according to 1)21(1 '). 

The refraction index (|25|l reaches its maximum for transversal propagation 
(fc|| =0,-d = tt/2) 

For _B ~ 10 • Bcr the declination of the refraction index (|26|l from unity 
exceeds that value for gases at atmospheric pressure in optic range by an 
order of magnitude, for B ~ 1000 • Bcr it reaches the value, characteristic of 
transparent liquids and glass; the refraction index H26|l becomes equal to that 
of diamond (n = 2.4) for B = 27 ■ lO^^^Gauss. 

Contrary to the just considered case of Mode 2, the polarization tensor 
eigenvalue K3 ((15|l does not contain a contribution linearly growing with the 
external field, as well as the resonance. For Mode 3 the dispersion equation 

with its r.-h. side given as 1(15(1 has the solution 

Po-H + kl^^, (27) 

where 

Z = l-ii(ln£-C-1.2l). (28) 

The known absence of asymptotic freedom in QED manifests itself in the 
negative sign in front of the logarithm in ((28() . This results in pathological 
consequences for the fields, as large as ~ Bcr exp(37r/Q!). In this domain the 
coefficient of fc^^ in 1(27(1 becomes, as the field grows, first less than zero and 
later greater than unity. The corresponding dispersion laws are nonphysical, 
since they lead to the group velocity greater than unity. In the negative slope 
case in 1(27(1 . e > i3exp(— 0.21 — C — 37r/a)/i?cr > 1, the dispersion curve 
enters the sector ((12(1 acceptable for the Cherenkov radiation. But this is the 
Cherenkov emission of tachyons! It is also odd, that in the latter case e.-m. 
waves may only propagate inside the cone < tan?9 < — 1 + a/SnZ with 
its axis along the external field, irrespective of the way they are produced. 
This domain of exponentially large external fields is not of our interest in the 
present paper. 
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For the fields that are not exponentially large, with the logarithmic terms 
of the order of unity, one should treat all the terms, marked by the coefficient 
a/37r in 1)27(1 . as small. Then, finally, the dispersion law for Mode 3 becomes: 

kl ^kl + kl{l^ ^) (29) 

Notably, the field-containing logarithmic terms have cancelled from here. 
Thus the dispersion law (|29|l of Mode 3 is saturated in the sense that, unlike 
Ea. ((21|l for Mode 2, it has reached the universal form, independent of the 
external field in the superSchwinger limit. The refraction index of Mode 3 
corresponding to is 

= 1 + sin^ ^. (30) 

DTT 

Again, the same as (|26f) . the maximum refraction in Mode 3 is achieved at 
perpendicular propagation, = 7r/2: 

= 1 + 3.8 • 10-"*. (31) 

This refraction index is of the order of that of gaseous ammonia and cannot 
be made larger by increasing the external field any further. 

4 Conclusion 

We have found that in the asymptotic case of external magnetic fields B that 
can be orders of magnitude larger than the Schwinger value of 4.4 x 10^'^Gauss, 
the refractive capacity of the magnetized vacuum grows unlimitedly with this 
field for electromagnetic radiation belonging to the polarization Mode 2, but 
comes to a saturation at a moderate level of corrections ~ a/Sir for Mode 3. 
For the "parallel energy" of the photon not close to the cyclotron resonance, 
/cq — A:| -C 4m^, the refraction effects for Mode 2 essentially exceed the above 
small corrections, typical for nonasymptotic domain, already for _B ~ 10 • Bcr. 
A regime is established in the range of photon frequencies/energies, extend- 
ing from zero and up to soft 7-rays, for which the dispersive properties of 
the magnetized vacuum are in each mode independent of the photon fre- 
quency/energy, but do depend upon the direction of its propagation. Apart 
from the fact that the refraction index in Mode 2 for the propagation non- 
parallel to the external field numerically grows with the field, it is remarkable 
that the angle between the group velocity and the direction of the photon 
momentum grows, too, the wave packet being attracted by the line of force 
of the external field. In this way the effect of 7-quantum capture by a strong 
magnetic field, known to exist due to resonance phenomena associated with 
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free and bound pair creation, is extended to lower energy ranges. Therefore, 

not only hard 7-rays, but also X-rays, light and radio-waves undergo strong 

dispersive influence of the magnetized vacuum, when the magnetic fields are 

of the order of magnitude of those estimated to exist in magnetars. In view 

of this fact the e-m energy canalization phenomena may become important 

not only within the traditional context described in Subsection 3.2 above, but 

also as applied to the scattering of e-m waves falling onto the magnetic field 


from outside . These may be, for instance, the X-rays emitted from the ac- 
cretion disk or from the pulsar surface outside the region where the magnetic 
field enters it. Recently the problem of the bending of electromagnetic radia- 
tion by the dipole magnetic field of a neutron star was addressed^", and the 
competition of this process with the effects of gravity was considered. (The 
author is indebted to Prof. H. Mosquera Cuesta who attracted his attention 
to that work in the course of Workshop). We insist, however, that such effects 
cannot be adequately treated disregarding the refraction index dependence on 
the direction of propagation and using the quadratic-in-the-field expressions 
for polarization operator, only valid in the low field limit, as is the case in 
Ref.^. 
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Appendix 

In this Appendix we derive t he asym ptotic expansion, presented in Subsection 
3.1, from expressions of 

The three eigenvalues Hi, i = 1,2,3 of the photon polarization operator 
in the one-loop approximation, calculated using exact propagator of electron 
in an external magnetic field, can be expressed as linear combinations of the 
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three functions Sj 

Kl = ~ {Zl + Zl) Si, K2 = (^1^2 + ZlY.\) , 

1^3 = (-22^3 + ^;iEi) , 

where the new notations for the momentum variables: 

Zl = k?i — fcg, Z2 = k\ 



(32) 



(33) 



are introduced, k"^ = Zi + Z2- Here Ej are dimensionless functions of the three 
ratios Bcr/B, Z2Bcr/Tn'^B, ziBcr/fn^B, equal to 



B 



sinht t->o sinht 



(34) 



(35) 



(2) / Bar Z2Bcr ZiBcr\ _ 2a 



B ' m^B ' m'^B 



„ , , d?7 - - ■ - X 
B I smhf 



where 



and 



exp -Z2 — -7; Zl , ^ t 1 - 1 



ai{t,r]) = 



eB 



cosh t — cosh tri 
2 sinh t ' 

1 — ?7 sinh(l + r])t 
~2 2 sinh i ' 



AeB 



1-rf 
<72{t,r]) = — - — cosht, 



cosh t — cosh r]t 
2sinh^t 



(36) 
(37) 

(38) 
(39) 
(40) 



newhepp: Workshop SMFNS/ICIMAF - February 1, 2008 



16 



The designation "lim" in (|35|l stands for the asymptotic hmit 

lim^ = i-^, z = l,2,3. (41) 
t^o sinht 4t ^ ^ 

The fact that do not depend on the fourth possible dimensionless variable 
z\l Z2 seems to be an approximation- independent manifestation of analyticity 
properties due to dispersion relations of Kramers-Kronig nature. 

Consider Y^i^ first. It does not depend on the photon energy and mo- 
mentum. With the use of the notation 

g,{t) = J ^a,{t,Tj) dr;, (42) 

can be represented as 

The integrals f I42|l are explicitly calculated to give 

1 / sinh2i \ , , cosht 



1 / , sinht 

TT- cosht 

sinh-^t V t 

Our goal now is to find the asymptotic behavior of (|43|l as 

B 

>oo. 

JDr/r 



53(t)-— 727 (cosh t-^^). (44) 



(45) 



The integrals in l|43(l do not contain singularities in the integrands at i = 0, but 
would diverge at t ^ oo if one just sets the hmiting value exp{—tBcr/B) — 1 
in the integrand. For this reason we should divide the integration domain into 
two parts . Besides, (8(72(^)7 sinht) — > 1 when t 00 and, hence, we have to 
add and subtract this limit beforehand in the integrand of ^2^^ Handling the 
cases of z = 1,3 does not require the like, since (51^3/ sinht) — > sufficiently 
fast. 
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B J \ smht J Jrp \ B J t 



T 



sinlit t / ./t \ sinht / B, 

+ ln(^^^ +C + lnT. (46) 

Here T is an arbitrary positive number, Si2 is the Kronecker delta and C 
is the Euler constant. We have omitted the exponentials from the first two 
integrals after the second equality sign in l|46|l . since now the resulting in- 
tegrals converge, and used the known asymptotic expansion of the standard 
exponential-inte gral function, which is - up to terms linearly decreasing with 
B/Bcr - equal tcP3 

/ tB„\ dt ( B„ \ , , 

, exp(-^j-^(ln^+lnT + cj (47) 

The most slowly decreasing term neglected from 
3/^dtexp(-tBcr/S)g^(t)/sinht is (5,;i(3B„/4B) ln(Bcr/S), since 
(7i(t)/sinht ~ (l/4t^) as t ^ oo. Other neglected terms decrease at 
least as fast as Bcr/B, since {'igi{t) / smht) — 5i2 decreases exponentially, like 
~ exp(— 2t), for i = 2, 3 when t is large. 

Numerical calculations using MATHCAD code give for the constants, 
involved in gHJ (d/i^/dT = 0) 

^ /3g.(t) i\ ^ r ng^{t) 



the following values: 

hi = 1.21, h2 = -0.69, hs = 0.21. (49) 
Finally, in the asymptotic regime B / B^r ^ 1, one has 
^(1) 2a ( , Brr „ , B 



S- = -^ln^ + C + ;.. + -fej (50) 

with the accuracy to terms, decreasing at least as fast as integral powers of 
the ratio Bcr/B and to the slower term {a / 2tt) {B / B)\n{Bcr/ B), omitted 
from y!^^\ 

Turn now to sf ^ Eq. This depends upon the three arguments, as 

indicated in H36|l . We shall be interested in the asymptotic domain, depicted 
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by the condition H45|) and 



2eB 

Zl 



(51) 



As for the ratio 2eB / Z2 we shall keep it finite, whenever it makes sense. 

The asymptotic expansion of (|38|) . 1)39(1 . 140(1 in powers of exp(— and 
exp(ryi) produces an expansion of ((^{6(1 into a sum of contributions coming from 
the thresholds ((16(1 . the singular behavior in the threshold points originat ing 
from the divergencies of the t-integration in l(3()() near t = oo (see l ^ l '^l l ^^ l for 
detail). The leading terms in the expansion of 1(38(1 . 1(39(1 . ((40(1 at i ^ oo are 



sinht 

0"2(t, -q) 
smht 



1-77 



exp(t(77-l))), 



(l + 2exp(-2t)), 



(52) 
(53) 



(73(^, ??) 

sinht 



2exp(-2i). 



(54) 



With the change of the variable t — t/eB and taking into account that it 
follows from ((37(1 that M (0x3,77) — 1/2, one obtains for 1(36(1 near the lowest 
singular thresholds {n — 0, n' — 1 or n' = Q, n — 1 for i = 1 in ((16(1 . ri = ri' = 1 
for 7 = 3, and n = 71' = and n = 77' = 1 for 7 = 2) 



^(2) ^ 2aeS 



^eB r\ ^ . 2 N 

drj dr exp(— 777 r) 

TT J-i Jo 



ai{eBT, ri) 



sinh ei?r 



exp 



Z2 21(1-772) 



2eS 



After the integration over t we get, e.g., for E]^' 



(2) 



.(2) 



AaeB 



Arj (1 — 77) X 



(55) 



exp I 



Z2 ' 

2eS / 



Itti^ + 4(1 - 77)eS + zi(l - 77^ 



47n2 



1 



4(1 - 77)eS 



(56) 

The pole in the above expression, caused by the integration over i, turns into 
the inverse square root singularity after the integration over 77 (c/. derivation 
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of ((T7|l from ifTHl l. In the limit (45), (51), when B » m^, i? > \zi\, no singu- 
larity remains in this expression (it is shifted to infinitely remote region) and 
we are left with S-j^ ' — (2a/7r)(exp(— Z2/2ei3) — 1). The same situation oc- 

(2) (2) 

curs for S3 , and for higher thresholds (also for contributions into S2 other, 
than those coming from the first term in H53|) 'l. The result of the calculation 
analogous to (|56|l for Eg^^ is (4a/7r)(exp(— 22/2eB) — 1). We conclude that in 
the limit (45), (51) there are no cyclotron resonances in the eigenvalues ki^ 3 
according to (|32|l . also that E^^'' does not introduce a singular contribution 
into K2- Consequently, there is no reason to keep the ratio eB / Z2 finite, when 
B ~> 00, in Ei^ 3, since Z2 may infinitely grow on the dispersion curve only 
when there is a resonance. 

Therefore, we should consider only the limit, when all the three arguments 
in Til 3 tend to zero. Handling this limit in H3t)|l is straightforward: 

The both integrations here converge, and hence this contribution decreases 
like zi/eB and Z2/eB when B 00. This is to be neglected within our scope 
of accuracy. 

(2) 

The situation is different with E2 . The resonance behavior is present 
here due to the contribution of the leading asymptotic term (1 — 77^)/4 in H53|l . 
It is responsible for the first threshold at —zi — 4to^ ( the ground Landau 
state n = n' = in H16|l ). which remains in its place when B ~* 00. Hence we 
should keep the ratio Z2/ eB nonzero when passing to the limit of large fields 
(since Z2 ^ 00 near the singular threshold on the dispersion curve) for the 

l2\ 

contribution of this term into Eg . The contributions of nonleading terms in 
the expansion 153|l into Ej^"* are nonsingular and should be treated along the 
same lines as El. 3 above. They decrease like zi/eB, Z2/eB to be neglected. 
Finally, for (|^ we are left in the limit (|^ with 



E(^)^^rdt expf-^^ t d.^x 



^ Jo \ B J 4 

e.J-Z2^-z,'-^t]-l 
^\ eB ' 4eS / 



(58) 



By the change of the integration variable t = eBr and the use of the asymp- 
totic form M{eBT,ri) = 1/2, eBr ^ 1 we finally obtain - after the r- 
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integration - the leading contribution into in the limit H45() , H51() 



By combining Eas. 1)59(1 and ((5()|l according to ((32|l . and bearing in mind 
((33|l and ((34|l , we come to the expressions , ((14|l , ((15|l for the polarization 
operator eigenvalues, given in the body of the paper above, provided that we 
also neglect the constant and logarithmic terms in K2 coming from Ef^ dSOll 
as compared to the terms growing linearly with B. 

References 

1. S. L. Adler, J. N. BahcaU, G. G. CaUan and M. N. Rosenbluth, Phys. 
Rev. Lett. 25, 1061 (1970). 

2. A. E. Shabad, Lettere al Nuovo Cimento 2, 457 (1972). 

3. A. E. Shabad, Ann. Phys. 90, 166 (1975). 

4. A. E. Shabad and V. V. Usov, Nature 295, 215 (1982), Astrophys. Space 
Sci. 102, 327 (1984). 

5. V. V. Usov and A. E. Shabad, Pis'ma v ZhETF 42, 17 (1985) {Sov. 
Phys.- JETP Letters 42, 19 (1985)). 

6. H. Herold, H. Ruder and G. Wunner, Phys. Rev. Lett. 54, 1453 (1985). 

7. A. E. Shabad and V. V. Usov, Astrophys. Space Sci. 117, 309 (1985), 
128, 377 (1986). 

8. V. B. Bhatia, Namrata Chopra and N. Panchapakesan, Astrophys. J. 
388, 131 (1992). 

9. V. V. Usov and D. B. Melrose, Aust. Journ. of Phys. 48, 571 (1995). 

10. V. V. Usov, ApJ. L87, 572 (2002). 

11. M. G. Baring, Astrophys. J L69, 440 (1995). 

12. A. C. Harding, M. G. Baring and P. L. Gonthier, Astrophys. J. 246, 
476 (1997). 

13. M. V. Chistyakov, A. V. Kuznetsov and N. V. Mikheev, Phys. 
Lett. B434, 67 (19 98); M. V. Chistyakov and N. V. Mikheev, 
arXiv:hep-ph/010172l| r. 

14. Andreas Reisenegger, The present volume (2003). 

15. V. M. Kaspi, Pulsar Astronomy - 2000 and Beyond. ASP Confer- 
ence Series 202 (1999) (M.Kramer, N.Wex and R.Wielebinski, eds.), 
arXiv:astro-ph/9912284 

16. V. V. Skobelev, Izv. Vyssh. Uchebn. Zaved. Fiz. NolO, 142 (1975); Yu. 
M. Loskutov and V. V. Skobelev, Phys. Lett. A56, 151 (1976). 



newhepp: Workshop SMFNS/ICIMAF - February 1, 2008 21 




(59) 



17. A. E. Shabad, Kratkie Soobtchenia po Fizike ( Sov. Phys. - Lebedev Inst. 
Reps.) 3, 13 (1976). 

18. I. A. Batalin and A.E.Shabad, ZhETF 60, 894 (1971) {Sov. Phys.- JETP 
33, 483 (1971)). 

19. H. Perez Rojas and A. E. Shabad, Ann. Phys. 121, 432 (1979). 

20. H. Perez Rojas, ZhETF 76, 3 (1979) (Sov. Phys.- JETP 49, 1 (1979)); 
H. Pe'rez Rojas and A. E. Shabad, Ann. Phys. 138, 1 (1982). 

21. A. E. Shabad, Trudy Fizicheskogo Instituta im. P.N. Lebedeva 192, pp.5 
- 152 (1988). 

22. A.E.Shabad, Polarization of the Vacuum and a Quantum Relativistic Gas 
in an External Field (Nova Science Pubhshers, New York, 1991) Proceed- 
ings of the Lebedev Physics Institute, Academy of Sciences of the USSR, 
Series editor N.G.Basov, Volume 191, edited by V.L.Ginzburg. 

23. V. I. Ritus, ZhETF m 1517, (1975); 73 807, (1977). 

24. Wu-yang Tsai, Phys. Rev. D. 10, 2699 (1974); V. N. Baier, V. M. Katkov 
and V. M. Strakhovenko, ZhETF 68 403, (1975). 

25. Guojun Qiao, The present volume (2003). 

26. V. L. Denisov, I. P. Denisova and S. I. Svertilov, Doklady Akademii Nauk 
380, 754 (2001); arXiv:astro-ph/0110705, 

27. Handbook of Mathematical Functions, edited by M. Abramovitz and I. 
Stegun, (National Bureau of Standards, New York, 1964). 



newhepp: Workshop SMFNS/ICIMAF - February 1, 2008 22 



